function sim = set_up_mit_entry(eq_SS, param, glob, options, eq_mit)

% entry shock
  %  load ../../initial_guess/initial_guess.mat
    
   sim = eq_mit;
   
   sim.D = eq_mit.D/eq_mit.D(end)*eq_SS.D;
   
    sim.A          = zeros(options.T, 1);
 %   sim.A(1)       = -.01;
 %   sim.A(2:100)   = sim.A(1)*(.5.^(2:100));
    sim.A          = exp(sim.A);
       
    sim.cE         = param.ce*ones(options.T,1);
    sim.cE(1)  = param.ce*(1+0.2);
  %  sim.cE(1:100)  = param.ce*(1+.1*(0.3).^((1:100)-1));
    
    sim.M          = param.M*ones(options.T,1);

    sim.cF         = param.mu_f*ones(options.T,1);

%     sim.C          = ones(options.T, 1);
%    sim.C(1:50)     = sim.C(1:50).*guess(:,1);
     
%    sim.W          = ones(options.T, 1);
%     sim.W(1:50)     = sim.W(1:50).*guess(:,3);

%    sim.D = eq_SS.D*ones(options.T,1);
%     sim.D(1:50)     = sim.D(1:50).*guess(:,2);

    sim.mu         = nan(size(eq_SS.L,1),options.T);
    sim.mu(:,1)    = eq_SS.L;
    
    sim.y          = nan(size(eq_SS.L,1),options.T);
    sim.p          = nan(size(eq_SS.L,1),options.T);
    sim.bp         = nan(size(eq_SS.L,1),options.T);
    
    sim.enter      = nan(size(glob.agrid,1), options.T);
    sim.exit       = nan(size(eq_SS.L,1),options.T);
    
end